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Abstract 



The thermal conductivity, k, of the heavy fermion superconductor UPts 

was measured down to Tc/10. The absence of a hnear term in the temperature 

dependence as T— >0 strongly suggests there are no zero energy quasiparticle 

excitations, in contradiction with the gapless behaviour often inferred from 

specific heat. A non-vanishing anisotropy ratio Kc/'^b as T^O establishes a 

new property of the gap structure: the presence of nodes along the c-axis. 

Furthermore, recent calculations by Fledderjohann and Hirschfeld show these 

cannot be point nodes with a linear k-dependence. 
PACS numbers: 74.70. Tx, 74.25.Fy 
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One of the central endeavours in the field of unconventional superconductivity is the 
precise determination of the gap structure. Numerous efforts are currently being made to 
establish whether or not a d-wave gap is realized in the cuprate superconductors. In the 
heavy fermion superconductor UPta, a number of candidate states are still possible after 
10 years of investigation. Any state with a line node in the basal plane of the hexagonal 
crystal structure, of which there are more than ten, is compatible with most of the existing 
data [1,2]. In particular, any such state is thought to be compatible with the compelling 
anisotropy of transverse sound attenuation [3]. A related issue of interest is the impact of 
impurity scattering on the low temperature properties of such superconducting states. The 
appearance of gapless behaviour, i.e. a residual density of quasiparticle excitations at T=0, 
even for small concentrations of non-magnetic impurities is under active investigation in the 
high Tc cuprates. In heavy fermions, no detailed study has been made of this issue - even 
though the specific heat of UPta, for example, has often been held as indicative of a residual 
normal fiuid in this compound - nor have there been quantitative comparisons with existing 
calculations. 

The detailed nodal structure of the gap and the effect of impurity scattering will bear 
directly on the various scenarios for the multicomponent phase diagram of UPts [2] . Both 
properties can only be studied at very low temperatures, typically in the region of Tc/10. 
In a previous paper [4], we showed how the thermal conductivity is a privileged probe of 
quasiparticles in UPts. In particular, it is sensitive to the anisotropy of the gap structure. 
Motivated by the recent calculations of Fledderjohann and Hirschfeld [5], we have now 
measured the thermal conductivity of UPta down to Tc/10 and report two new facts: 1) 
there is no trace of a residual normal fiuid at T=0 (at the 1% level) and 2) the gap must go 
to zero along the c-axis, in addition to the well-established line node in the basal plane. A 
comparison with the results of Fledderjohann and Hirschfeld further suggests that the gap 
cannot vanish linearly at a point, thereby weighing against the Eig representation for the 
superconducting order parameter of UPta [2] . 

The details of the experiment and the crystal characteristics are described in Ref. [4]. 
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The thermal conductivity k{T) of UPts was measured for a heat current along each of 
the two high-symmetry directions of the hexagonal lattice (J//b-axis and J//c-axis), as a 
function of temperature down to Tc/10~50 mK. The low-temperature results are shown in 
Fig. 1, plotted as k/T vs T, while the overall behaviour is shown in the inset. The normal 
state behaviour below T+=0.5 K was obtained by applying a magnetic field of 3 T > Hc2(0) 
(dashed lines). In this paper, we will only be concerned with the low temperature- low field 
phase B, which exists below T7=0.44 K [2]. 

In order to use heat conduction as a measure of electronic transport one must ensure 
that phonons do not contribute significantly, i.e. that Kph«Ke in the total conductivity 
K—Kph+Kg. The safest estimate of the maximum possible phonon contribution is obtained by 
using the formula Kph—Cphyph-^ph/'^, where Cph—PT^ is the low temperature phonon specific 
heat, Yph is the average sound velocity in the direction of the heat current and Aph is the 
phonon mean free path. One then assumes that Ap^ takes on its maximum value, namely 
the size of the crystal, equal to 0.7 mm [4]. Prom pubhshed data, P—20 J m"^ [6], in 
agreement with the known sound velocities [6]. The correct average velocity Vp/j [7] is 1880 m 

for acoustic waves traveUing along the b-axis and 1440 for the c-axis [6]. The maximum 
heat conduction by phonons in our crystal is therefore Kph—S5 (67) T^ mW K~^cm~^ for the 
b (c) axis. Por T<150 mK, this represents at most 15% of the measured Hh and 6% of Kc- 
Of course, scattering by quasiparticles, which are present in sufficient numbers to account 
for more than 85% of will decrease Ap^ from this maximum possible value. This decrease 
will occur very rapidly as the temperature is increased, with Aph reaching a mere 5 fim at 
T=Tc [4]. Therefore, Uph is certainly less than the upper bounds of 15% and 6%, and can 
thus be neglected. We point out that while in our crystal the electronic mean free path is 
long enough to ensure that K,ph«Ke, this is less likely to have been the case in previous 
measurements on polycrystalline UPts, where Ke was 4 to 6 times smaller [1], not to mention 
measurements on other heavy fermion compounds [1]. As a result, we believe this is the first 
time heat transport by heavy fermion quasiparticles is reliably measured down to Tc/10. 
This allows us to examine the possibility of a residual normal fluid at T=0 and it provides 
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us with a new and powerful probe of the gap structure. 

Ideally, the question of a residual normal fluid should be answered by low-temperature 
measurements of the specific heat. Unfortunately, this has proven difficult both in UPta and 
in the high Tc cuprates, for two reasons. First, a number of sizeable non-linear contributions 
to C(T) exist at low temperature. In UPta, the quasiparticle contribution, although large, 
is completely overwhelmed below 100 mK by a huge upturn in C/T of ill-understood origin 
[8] . Secondly, even if a residual linear term is extracted reliably, it cannot automatically be 
attributed entirely to fermion excitations, for other mechanisms also lead to a linear specific 
heat. For UPts, the standard approach has been to extrapolate down from the roughly 
hnear behaviour of C/T observed above about 100 or 150 mK [1,2,8], assuming it to persist 
down to T=0. This procedure yields an intercept at T=0, called 70, which typically ranges 
from 10 to 40% of the normal state value ■jn in high quality samples [1,8]. In our crystal, 
70=16% 7iv [9]. 

In recent years, several theories have been proposed for this extrapolated 70, taken to be 
an intrinsic property of the quasiparticle fluid in UPta [1,10,11]. In light of the strong sample 
dependence, we are more inclined to regard this as an extrinsic property. For unambiguous 
evidence about a possible residual normal fluid at T=0, we therefore argue that one must turn 
to a more reliable measurement, such as heat conduction, free from the problems mentioned 
above. The observation of a finite k/T at T— >^0, such as seen in the high Tc compound 
YBa2Cu307, would be a direct indication of zero-energy quasiparticle excitations. It is clear 
from Fig.l that a smooth extension of the k/T data to T=0 leads to a negligible intercept. 
Of course, strictly speaking, an abrupt flattening off of both curves below 40 mK cannot 
be excluded; in this case the largest linear term would represent 2% (4%) of its normal 
state value at T=0, equal to Lo/po [4], for J//c (J//b). We therefore conclude that, in the 
absence of any real sign of a residual linear term in our data, there is no evidence for gapless 
excitations in UPta. 

More quantitatively, we note that at our lowest temperature T=0.1 T^, k/kn=3% (on 
average), whereas an extrapolated C/T for the same crystal gives 7(T=0.1Tc)/7Ar=30%. 
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Could it be that the mean free path due to impurity scattering in the superconducting 
state is 10 times smaller than in the normal state? One of the few measurements of the 
quasiparticle mean free path below Tc comes from the degradation of a heat current by 
vortices. The increase in the thermal resistivity 1/k with increasing magnetic field just 
above Hd is proportional to the density of vortices B/$o; the carrier mean free path Iq and 
the vortex effective diameter D. From their measurement of this effect on a crystal of a 
quality similar to ours, Behnia et al. [12] obtained IqD ~ 3000 nm^ at T=0.2 K. If D is 
taken to be the coherence length, then D ~ ^(T = 0.2K) ~ — 10 nm [1,6], and Iq ~ 300 
nm. As we will see below, this is very close to the mean free path we estimate for the normal 
state. Therefore, it appears unlikely that the relaxation time in the superconducting state 
is suppressed by an order of magnitude. 

The whole question of the impact of impurity scattering on unconventional gap struc- 
tures, and vice versa, was treated by several authors in the mid-80s (see Refs.[l,13,14] and 
references therein). Within a weak-coupling ECS theory, Hirschfeld et al. [13] showed that 
a self-consistent treatment of impurity scattering for a gap with line nodes can lead to a 
residual density of quasiparticle states, showing up as a finite C/T and k/T at T=0. Re- 
cent calculations of this kind were performed by Fledderjohann and Hirschfeld [5] for three 
uniaxial gap structures, each with a line node in the basal plane: 1) a polar gap (with no 
other nodes) and 2) two hybrid gaps (with in addition a point node at each pole, i.e. along 
the c-axis). One of the hybrid gaps, which we call hybrid-I, vanishes linearly in k at the 
point nodes, while the other, called hybrid-ll, vanishes quadratically. Some of the states 
allowed by hexagonal symmetry to which these gaps correspond are listed in Table I. The 
calculations so far assume a single ellipsoidal Fermi surface and s-wave scattering, and they 
require two input parameters, the impurity scattering rate Fq and the scattering phase shift 
Sq [5] . Much information about the nature of the scattering can be obtained from the normal 
state properties. From the fact that the anisotropy of transport (both heat and charge) is 
independent of temperature, i.e. that Kc/Kiy=2.8 and (Jc/cTfe=2.7 all the way from T=0.1 K 
(in the normal state) where elastic scattering dominates to T=0.8 K where inelastic electron- 
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electron scattering dominates [4], it is very likely that the anisotropy of 2.7-2.8 is due to the 
Fermi velocity, and both impurity-electron collisions and electron-electron collisions should 
be well-described by isotropic (s-wave) scattering [5]. From de Haas- van Alphen measure- 
ments (dHvA) [15], the Fermi surface is known to be made of several sheets, and a single 
ellipsoid is certainly an oversimplification. However, the usual Dingle plot analysis yields 
fairly uniform scattering rates, with r = 1/TdHvA = 2 - 4 x 10~^^ sec, in crystals of a quality 
comparable to ours [15]. In temperature, this corresponds to TdHvA — 0.2 - 0.4 Tc. The scat- 
tering rate Fq appropriate for transport, which is less sensitive to small-angle scattering and 
dephasing, will be smaller than TdHvA by a factor which depends on the type of scattering. 
A value of Fo=0.1 Tc seems quite reasonable, corresponding to Iq — vf/Tq — 400 nm [15]. 
A separate estimate, obtained from the shear viscosity at Tc [14], gives F(Tc)=ro(p(Tc)/po) 
= 0.22 Tc, so that again Fq — 0.1 Tc. Finally, Fledderjohann and Hirschfeld took 5o=7t/2, 
the unitary limit, seeing as weaker scattering leads to sharp disagreement with observed 
properties [5]. Theoretical arguments for such a limit can be found in Refs. [13] and [14], 
and in references therein. 

The results of the calculations for Hf, [5], for two values of Fq (0.1 and 0.01 Tc) are 
reproduced in Fig. 2, alongside our data, as a function of reduced temperature T/Tc. The 
data are plotted as fi;/fi;Ar=fi;(H=0)/fi;(H>Hc2), where fi;Ar=T/(a-|-bT^) [4] and Tc=Tc =0.44 
K. Given that b/a=0.25 [4], deviates from linear behaviour (elastic scattering) by 
only 7% at T=0.3Tc. This means that inelastic scattering can safely be neglected in the 
calculations below that temperature, shown therefore as k/T normalized to 1 at Tc. Only 
the curves for a polar and a hybrid-II gap structures are shown; the corresponding curves 
for the hybrid-I gap lie in between [5]. Inspection of both data and calculated curves reveals 
that the rapid increase in k/kat with temperature at such low T/Tc - an order of magnitude 
faster than in a conventional superconductor [4] - is well reproduced by the calculations; 
this is a convincing confirmation of a line node in the basal plane of the gap structure 
of UPt3 (phase B). However, the data show less curvature than cither of the curves with 
Fo/Tc=0.1 and, indeed, will not smoothly extrapolate to any significant intercept at T=0, 

6 



such as expected from the theory. In this sense, the observed behaviour is more compatible 
with calculations based on a much smaller Fq, say 0.01 Tc, with a temperature dependence 
close to that of the polar gap, although with a magnitude closer to the hybrid-II gap. Until 
the real Fermi surface is used both in the calculations and in the estimates of a transport 
Fq, related directly to the measured po, it is difficult to make firm conclusions from this 
comparison. Nevertheless, it does seem as though an unreasonably low scattering rate is 
needed to keep the number of zero energy quasiparticle excitations obtained in the current 
self-consistent calculations at the low level observed in the experiment. 

The anisotropy of electronic heat conduction has long been known to be a useful probe 
of gap anisotropy, even in s-wave superconductors. In Ref.[4] we showed how the ratio Kc/ 
is a direct probe of the anisotropy of the gap in UPts, insofar as it is constant above Tc and 
starts falling immediately below Tc (see inset of Fig. 3). The fact that it decreases rather 
than increases suggests there are more thermally excited quasiparticles with velocities along 
the b-axis than along the c-axis. This could either result from a finite gap being larger 
along c than along b (anisotropic s-wave gap), from the presence of nodes in the gap along 
the b-axis in the absence of any along the c-axis (polar gap), or from the presence of nodes 
along both axes (hybrid gap) provided the nodal structure is such that more quasiparticles 
have v/ /b. This, for example, would exclude an axial gap (with only point nodes along 
the c-axis) [4]. Note, however, that the current discussion relies on the assumption that 
the scattering rate does not change with temperature in some unexpected way below Tc. 
To be free from such ambiguity, the analysis must be done in a regime where the strong 
electron-electron scattering is not important, namely below 150 mK or so. Moreover, it is 
in the limit of T— >0 that a measurement of Kc/ Kf, becomes particularly useful. Indeed at 
T<<Tc, the regions of the gap very close to the nodes dominate the thermal properties and 
a measurement of heat conduction can then shed light on the detailed structure of the gap 
in the vicinity of both high-symmetry directions. These are the main justifications for the 
present low temperature study. 

Our results for Kc/ kj, below 0.3 Tc are shown in Fig.3, normalized at Tc=T~ (i.e. divided 

7 



by 2.8, the constant value of the normal state anisotropy). The striking finding is that the 
ratio does not go to zero as T^O. This definitively excludes not only an axial gap, but also 
a polar gap. Indeed, the latter will give Kc/nb—^O as T— >0, as a result of the clear difference 
between excitation of quasiparticles with v//c (across a finite gap) and with v//b (in the 
vicinity of a fine node) . We conclude that the gap of UPta-phase B must have nodes along 
the c-axis. No previous experiment imposed such a requirement on the gap structure of 
this material, except perhaps the V-shaped gap features observed in recent studies of point- 
contact spectroscopy [16]. Indeed, all other measured properties are so far compatible with 
a polar gap [1]. 

A comparison with the calculations of Fledderjohann and Hirschfeld [5] allows us to go 
further and gain insight into the specific nodal structure near k^—ky—O. Their results on 
Kc/ i^b for the three uniaxial gaps with line nodes (polar, hybrid-I and hybrid-II) are shown in 
Fig. 3. In the absence of gapless behaviour (i.e. for Fq = 0.01 Tc), they find that Kc/ k^-^Q as 
T— >^0 not only for the polar gap, as expected, but also for the hybrid-I gap. In essence, a gap 
vanishing at a point node with linear k-dependence does not cause as many quasiparticles to 
be excited thermally as a gap vanishing along a line. Remarkably, these authors found this 
not to be true for a point node with quadratic k-dependence, and the hybrid-II gap yields 
a finite value for Kc/ Hb as T^O. More specifically, such a gap on a spherical or ellipsoidal 
Fermi surface leads to no change in the anisotropy of heat conduction with respect to the 
normal state [5]. Note that the perfect isotropy will not be preserved for the real Fermi 
surface [5]. 

We conclude that the gap structure of phase B is unlikely to be of the hybrid-I type, as 
has been widely believed over the past few years [2]. Instead, our low temperature results 
favour a gap of the hybrid-II type, being the only uniaxial gap with the correct limiting 
behaviour of Kc/ Kb (see Table I). Of course, states with non- uniaxial symmetry should also 
be considered in future calculations. In the context of the 2D theory for the multicomponent 
phase diagram of UPts (see Refs. [2] and [17], and references therein), this appears to 
disqualify the (l,i) state of the Eig, Fi2g and Ei„ representations as possible candidates for 
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phase B, and leave only the (l,i) state of the representation with d-vector parallel to 
the c-axis [17]. 

In summary, measurements of the thermal conductivity of UPta down to Tc/10 have 
shed light on two important aspects of unconventional superconductivity: the possibility 
of a gapless behaviour and the nodal structure of the gap function. In the absence of any 
detectable linear term in k as T^O, we find no evidence for an intrinsic gapless regime, 
such as previously inferred from the specific heat. A comparison with calculations based 
on resonant impurity scattering and reasonable assumptions for the scattering strength and 
anisotropy leads to a significant discrepancy in the limiting value of Kb/T. The unusual 
observation of a finite value for the anisotropy ratio Kc/ as T— >-0 leads to new information 
about the gap structure of UPta in phase B: the gap vanishes along the c-axis, and it does so 
with a special k-dependence, not compatible with a linear point node but probably so with 
a quadratic point node. This places the most stringent constraints so far on the possible 
superconducting states, and hence on the various theoretical models, for UPta. 

We are grateful to A. Fledderjohann and P.J. Hirschfeld for extensive discussions and 
for allowing us to reproduce their calculations here. This work was funded by NSERC of 
Canada and FCAR of Quebec. L.T. acknowledges the support of the Canadian Institute for 
Advanced Research and the A. P. Sloan Foundation. 
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FIGURES 

FIG. 1. Low-temperature thermal conductivity of UPta, divided by temperature, for a heat 
current along the c-axis (open circles) and the b-axis (solid circles). Inset: k/T up to 0.8 K. The 
normal state behaviour, obtained by applying a field above Hc2, is also shown (dashed lines). 

FIG. 2. Thermal conductivity along the b-axis below 0.3 Tc, where Tc=T~=0.44 K. Upper 
frame: data plotted as k/kn, where «;Ar=T(a-|-bT^)~^ (see text). Lower frame: k/T calculated 
for two uniaxial gaps with line nodes (polar and hybrid-II) and for two values of the impurity 
scattering rate Fq, 0.1 Tc (solid lines) and 0.01 Tc (dashed lines) (after Ref.[5]). 

FIG. 3. Low-temperature behaviour of the anisotropy ratio Kc/k^, plotted as a function of 
reduced temperature and normalized to unity at Tc. The data (points) are compared with calcula- 
tions for three gap structures with a line node in the basal plane, with ro=0.01 Tc (after Ref.[5]). 
The dashed line is a linear fit to the data below 0.3 Tc. Inset: unnormalized data up to 0.8 K. 
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TABLES 

TABLE I. The gap structure of uniaxial states allowed by hexagonal symmetry (for strong 
spin-orbit coupling), and the limiting value of the anisotropy ratio Kc/ Kb expected as T— >0, in the 
absence of gapless behaviour (see Ref. [5]). The nodal structures include a gap going to zero at a 
point along the c-axis, either with a linear (LP) or a quadratic (QP) k-dependence, and along a 
line in the basal plane (line node). Only odd-parity states with d//c are listed (see Ref. [17]). 



Gap Nodal structure Superconducting states nd nb as T— > 

Anisotropic s-wave none Aig (A^ < Ac) 

Axial I LP Eiu (l,i) oo 

Axial II QP E2g{l,i) oo 

Polar line node Ai^ 

Hybrid I line node LP ^ig(l,i) 

Hybrid II line node QP 1 
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